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ABSTRACT 

We present predictions for number counts and redshift distributions of galaxies de¬ 
tectable in continuum and in emission lines with the Mid-infrared (MIR) Instrument 
(SMI) proposed for the Space Infrared Telescope for Cosmology and Astrophysics 
(SPICA). We have considered 24 MIR fine-structure lines, four Polycyclic Aromatic 
Hydrocarbon (PAH) bands (at 6.2, 7.7, 8.6 and 11.3/jm) and two silicate bands (in 
emission and in absorption) at 9.7/im and 18.0/im. Six of these lines are primarily 
associated with Active Galactic Nuclei (AGNs), the others with star formation. A 
survey with the SMI spectrometers of 1 hour integration per field-of-view (FoV) over 
an area of 1 deg 2 will yield 5 a detections of ~ 140 AGN lines and of ~ 5.2 x 10 4 
star-forming galaxies, ~ 1.6 x 10 4 of which will be detected in at least two lines. The 
combination of a shallow (20.0deg 2 , 1.4 x 10 _1 h integration per FoV) and a deep 
survey (6.9 x 10 -3 deg , 635 h integration time), with the SMI camera, for a total of 
~1000h, will accurately determine the MIR number counts of galaxies and of AGNs 
over five orders of magnitude in flux density, reaching values more than one order 
of magnitude fainter than the deepest Spitzer 24 /tm surveys. This will allow us to 
determine the cosmic star formation rate (SFR) function down to SFRs more than 
100 times fainter than reached by the Herschel Observatory. 

Key words: galaxies: luminosity function - galaxies: evolution - galaxies: active - 
galaxies: starburst - infrared: galaxies 


1 INTRODUCTION 

Studying the co-evolution of star-formation and black-hole 
accretion is one of main scientific goals of the SPace In- 
fraRed telescope for Cosmology and Astrophysics (SPICAR. 
SPICA will be equipped with two main instruments: th e 
SpicA FAR infrared Instrument (SAFARI; iRoelfsema et al.l 


* matteo.bonato@oapd.inaf.it 

1 http://www. ir.isas.jaxa.jp/SPICA/SPICA_HP/index-en.html 


2012h and th e SPICA Mid-infrared Instrument (SMlR. In 
Bonato et al.l (l2014al lbh we have presented detailed predic¬ 
tions for the number counts and the redshift distributions 
of galaxies detectable in blind spectroscopic surveys with 
SAFARI, accounting for both the starburst and the Active 

2 https://home.sron.nl/files/LEA/SAFARI/spica_workshop_2014 
/SMI_factsheet2.pdf 

https://home.sron.nl/files/LEA/S AFARI/spica_workshop_2014 
/KanedaH_SPICA_workshop_2014.pdf 
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Figure 1. Euclidean normalized differential number counts (on the left) and integral counts (on the right) at 24 jim given by the adopted 
model for galaxies as a whole (starburst plus AGN components; solid blue lines) and for AGNs alone (dotted violet lines). The vertical 
solid black and dashed cyan lines correspond to t he SMI camera detec ti on limits for the propo s ed deep and shallow sur v eys, respectively. 
Count s of galaxies (filled g reen circles) are from iTakagi et alj (12012 ) , I Clements et all (120 111) , iBethermin et all d201C)h , iLe Floc’h et al.l 
J 2 OO 9 I) . [shupe~ et alJ (|2008h , IPapovich et al.1 j2004f ) and lSanders et all (|2003h . The latter data show an excess at the highest flux den sities 
due to the enhance d number of bright galaxies in the Virgo supercluster. AGN counts (red triangles) are from pTreister et al.l (120061) and 
I Brown et al . ( 2006}) . 


Galactic Nucleus (AGN) components. Here we focus on the 
SMI. 

The SMI has two basic observing modes: the wide- 
held imaging camera mode and th e spectrometer mode with 
two detectors (Spec-S and Spec-L; IKataza et al.ll2012T ). The 
technical specifications for the SMI used in this work are: 
R=1000 spectrometers, FoV = 150” x3", wavelength ranges 
20 — 27 fj,m (Spec-S) and 27 — 37 /im (Spec-L); R=20 wide 
held camera, FoV = 5 ( x 5) wavelength range 20 — 37 fim. 
The spatial resolution (FWHM) varies from 1.4" at 20 /im 
to 2.6" at 37/.an. The line detection limit (lhr, 5 ct) is in 
the range 6 — 23 x 10 _2 ° W m -2 for the camera, 2 — 7 x 
ir 20 Wm- 2 for the Spec-S and 3 - 10 x ir 20 Wm“ 2 for 
the Spec-L. The point source continuum sensitivity (lhr, 
5 a) in a low background region increases from ~ 10 /t Jy 
at 20 /im to ~ 60 /iJy at 37 /im; at 24 /tm it is ~ 10.5 /tJy. 
The survey speed for the 5 a detection of a point source 
with a continuum hux density of 40 /tjy with the cam¬ 
era is ~ 7 arcmin 2 /hr; for the detection of the line hux 
of 2 x 10~ 19 Wm _2 it is ~ 4arcmin 2 /h for Spec-S and 
~ 2 arcmin 2 /h for Spec-L. 

The SMI instrument is crucial to enhance the outcomes 
of the spectroscopic surveys carried out with SAFARI. The 
R = 1000 spectrometer is needed to observe One struc¬ 
ture lines with a resolution similar to SAFARI, whereas 
the R = 20 wide held camera is essential to uncover star¬ 
forming galaxies in the four broad and very bright Poly¬ 
cyclic Aromatic Hydrocarbon (PAH) bands at 6.2, 7.7, 8.6 
and 11.3pm. The lines that will be detected can come either 
from star forming regions or from nuclear activity or from 
both. 

The SMI spectroscopy will allow us to exploit the rich 
suite of MIR diagnostic lines to trace the star formation 


and the accretion onto the super-massive black holes up to 
high redshifts through both blind spectroscopic surveys and 
pointed observations. The MIR lines detectable by the SMI 
provide excellent diagnostics of the gas density and of the 
hardness of the exciting radiation held. The ratios of two 
lines having similar critical density and different ionization 
potential allow us to estimate the ionization of the gas, while 
the ratios of two lines with different critical density and sim¬ 
ilar ionization potential provide estimate s of the gas density 
in the region (ISpinoglio fe Malkanlll992i ). A comprehensive 
discus sion of infrare d density ind i cators is given bv iRubinl 
d 19891) . As shown in ISturm et al.l (l2002l f , mid-IR line ratio 
diagrams can be used to identify composite sources and to 
distinguish between emission from star forming regions and 
emission excited by nuclear activity. These diagnostic dia¬ 
grams are constructed by plotting pairs of line ratios against 
each other (for example [NeVI]7.63/[OIV]25.89 and [NeVI]- 
7.63/[NeII]12.81), in which different types of regions can be 
easily separated and distinguished. Similar diagnostic dia- 
gr ams with different set s of we aker lines h ave b een proposed 
bv ISpinoglio fc Malkanl (Il992li and IVoifl (1 19921) . The multi¬ 
plicity of possible combinations of lines allows us to adapt 
these diagnostic tools to different redshift ranges. In this re¬ 
spect, the complementary wavelength coverages of SMI and 
SAFARI substantially enhances the potential of the SPICA 
mission. 

The SMI camera will also substantially improve our 
knowledge of source counts in the mid-infrared (MIR) region 
by extending them to much fainter flux density levels than 
achieved by Spitzer and reaching a much better statistics. 
This allows a considerable improvement of our understand¬ 
ing of the cosmic star formation history. 

In this paper we use the ICai et ahl (l2013i 'l evolution- 
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Figure 2. Number of predicted detections of galaxies as a whole, 
including the AGN component (upper blue points) and AGNs 
alone (lower violet points) within A log S = 0.1 bins (A log S = 
0.2 bins for faint AGNs) for the proposed shallow and deep sur¬ 
veys with the SMI camera. The error bars are the sum in quadra¬ 
ture of Poisson errors and of the variance due to clustering (sam¬ 
pling variance; see text). The latter is relevant only for the deep 
survey. The vertical solid black and dashed cyan lines correspond 
to the survey 5 a detection limits. 


ary model as upgraded bv lBonato et al.l (l 2014b ). The model 
deals in a self consistent way with the emission of galaxies 
as a whole, including both the starburst and the AGN com¬ 
ponent, and was successfully tested against a large amount 
of observational data. 

The plan of the paper is the following. In Section[2]we 
briefly summarize the adopted model for the evolution with 
cosmic time of the IR (8-1000 /im) luminosity function. In 
Section[3]we discuss imaging observations with the wide field 
SMI camera. In SectionUwe present the relations between 
line and continuum luminosity for the main MIR lines. In 
Section[S] we work out our predictions for line luminosity 
functions, number counts and redshift distributions within 
the SMI wavelength coverage. In Section[6]we discuss possi¬ 
ble SMI observational strategies. Section]?] contains a sum¬ 
mary of our main conclusions. 

We adopt a flat ACDM cosmology with matter den¬ 
sity n m = 0.32, dark energy density 12a = 0.68 and 
Hubble constant h = i/o/100kms^ 1 Mpc -1 = 0.67 

dPlanck Collaboration XVj|2014h . 


2 EVOLUTION OF THE IR LUMINOSITY 
FUNCTIONS 

The ICai et all (120131 ) model, adopted here, is based on a 
comprehensive “hybrid” approach that combines a physical 
model for the progenitors of the spheroidal components of 
galaxies (early-type’s and massive bulges of late-type’s) with 
a phenomenological one for disk components. The evolution 
of the former population is describ e d by an updated ver¬ 
sion of the physical model bv lGranato et al.l (|2004l . see also 
Lapi et al. 2006, 2011). In the local universe these objects 
are composed of relatively old stellar populations with mass- 
weighted ages ^ 8-9 Gyr, corresponding to formation red- 
shifts a iZ, 1-1.5, while the disc components of spirals and the 
irregular galaxies are c haracterized by signifi cantly younger 
stellar populations (cf. iBernardi et al.ll20ld . their Fig. 10). 


Thus the progenitors of spheroidal galactic components, re¬ 
ferred to as proto-spheroidal galaxies or “proto-spheroids”, 
are the dominant star forming population at z ^ 1.5, while 
IR galaxies at 2 1.5 are mostly galaxy disks. 


In the case of proto-spheroids, the model describes the 
co-evolution of the stellar and of the AGN component, al¬ 
lowing us to deal straightforwardly with objects as a whole. 
This does not happen for disk components of late-type 
galaxies, whose evolution is described by a phenomeno¬ 
logical, parametric model, distinguishing between the two 
sub-populations of “cold” (normal) and “warm” (starburst) 
gal axies. AGNs are treat ed as a separate population. Follow¬ 
ing (Bonato^et^L] ([2014b) , we have asso ciated them to the 
late-type galaxy populations using the IChen et al.l (120131) 
correlation between star formation rate (SFR) and black 
hole accretion rate. Both type 1 and type 2 AGNs, with 
rel ative a bundances, as a f u nction o f lum inosity, derived by 
lHasingeil (120081) (see lBianchi, Maiolino, fe Risalitill2012l . for 
a review) are taken i n to acc ount. 

The IChen et al.l (I 2 OI 3 I) correlation does not apply to 
bright optically selected QSOs, which have high accre¬ 
tion rates but are hosted by galaxi es w ith _SFRs ranging 
from very low to moderate. As in iBonato et al.l (12014b ) 
we reckon with th ese objects a d opting the best fit evolu¬ 
tionary model by ]Croom_gt_aL| (|2009) up to 2 = 2. As 
shown bv IBonato et al. ( 2014bll . this approach reproduces 
the observationally determined bolometric luminosity func¬ 
tions of AGNs at different redsliifts. At higher redshifts 
optic al AGNs are already accounted for by the ICai et all 
(l2013h model which also accounts for redshift-dependent 
AGN bolometric luminosity functions. 


3 SURVEYS WITH THE WIDE FIELD 

CAMERA 

The SMI camera can substantially extend the mid-IR flux 
density range observed so far (by the Spitzer satellite). Us¬ 
ing the counts yielded by the model, that are strongly con¬ 
strained by observational data (see the left panel of Fig. [T]), 
we find a 24 pm 5 a confusion limit of 4.8 x 10~ 5 mjy. For 
comparison, the deepest Spitzer surveys reached 30 pJy us¬ 
ing an ext racti o n tec hniq ue ba sed on prior source positions 
at 3.6 pm (iMagnelli et al.H~201ll) : the completeness of the re¬ 
sulting catalog, estimated via Monte Carlo simulations to be 
80%, is however difficul t to assess. As illustrated by Fig. 5 of 
iBethermin et al] (l20ld ) at fluxes lower than 8.0 x 10 -2 mjy 
the 24 pm counts are endowed with substantial uncertainties 
due either to poor statistics or to substantial corrections for 
incompleteness. 

With the current technical specifications (see Sect. [T]), 
reaching a 5 <7 limit of 3.0 /iJy, i.e. going at least one order 
of magnitude deeper than Spitzer, requires ~ 635 h per FoV. 
Observations of a single FoV (6.9 x 10 -3 deg 2 ) will be enough 
to detect hundreds of sources per 0.1 dex in flux density. A 
survey at this limit would resolve ~ 92% of the background 
estimated using the adopted model. A determination of the 
source counts over a very broad flux density interval can be 
achieved adding a shallow survey with a detection limit of 
0.2 mjy, an integration time of ~ 1.4x 10 -1 h per FoV, cover¬ 
ing an area of 20.0deg 2 . The global observing time (without 
overheads) is then 1046 h. 
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Figure 3. Integral counts in 4 PAH lines of star-forming galaxies over the SMI wide field camera wavelength range. The solid black and 
the dashed cyan vertical lines correspond to the SMI camera detection limits for the proposed deep and shallow surveys, respectively. 
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Figure 4. Redshift distributions of galaxies detected in at least 1 to 4 PAH lines by the proposed shallow (left panel) and deep (right 
panel) survey with the SMI wide-held camera. 


Not surprisingly, the deep survey would be very costly of 
time. Is there a sufficient scientific motivation for it, given 
that the 24 /im counts are found to be rapidly converging 
already at much brighter flux density levels? The simplest 
answer is that fully exploiting the potential of a new instru¬ 
ment that provides a large boost in sensitivity is a must since 
it is the most direct way for exploring the unknown, look¬ 
ing for the unexpected. On top of that, in this case going 
deeper guarantees important information on a still poorly 
understood aspect of galaxy evolution: the effect of feed¬ 
back on the star formation history of low mass galaxies at 
high redshifts, hence also on the build up of larger galaxies 
via mergers. 

It is generally agreed that the flatter slope of the faint 
end of the galaxy luminosity function compared to that of 
the halo mass function is due fe edback, but how the fe edback 
operates is still unclear (e.g.. Isilk fc Mamonl l2012h . Vari¬ 
ous flavors of feedback have been advocated, including re¬ 
ionization, supernova explosions, tidal stripping and more, 
but the role of each effect is not well understood. The high 
surface density reached by the deep survey (right-hand panel 
of Fig. [I]) implies that it will detect low mass galaxies, with 
low SFRs at substantial redshifts, as discussed at the end 
of this section. Thus it will provide direct information on a 
key aspect of galaxy evolution. 

The right-hand panel of Fig. [T| shows the predicted in¬ 
tegral counts at 24 /im of galaxies as a whole (starburst plus 
AGN component; solid blue lines) and of the AGNs alone 


(dotted violet lines). The proposed surveys are expected to 
detect ~ 1.1 x 10 s galaxies and ~ 1.0 x 10 4 AGNs (shallow 
survey), and ~ 1.3 x 10 3 galaxies and ~ 94 AGNs (deep sur¬ 
vey). Interestingly, the surface density of AGNs detected by 
the deep survey is essentially the sam e as that o f the deep¬ 
est C handra survey (4 million seconds; iBrandt fc Alexander! 
120151) in X-rays. The mid-IR AGN counts are crucial, among 
other things, to assess the abundance of heavily absorbed 
AGNs, missed by X-ray surveys, that can contribute an im¬ 
portant fraction of the high-energy X-ray background. 

Figure [2] shows the number of predicted detections of 
galaxies and AGNs in A log S' = 0.1 bins (AlogS = 0.2 
for AGNs in the case of the deep survey). Below 0.2mjy, 
where observational data are largely absent (see Fig. [T]), the 
proposed survey will detect from ~ 13 to ~ 177 galaxies per 
A log S = 0.1 bin, and from ~ 2 to ~ 23 AGNs per A log S = 
0.2 bin. At lOmJy, where the Spitzer counts have a very 
poor statistics (Fig. [T]) , the proposed survey will detect ~ 76 
galaxies and ~ 24 AGNs per AlogS = 0.1 bin. 

The error bars plotted in Fig.[2]include both the Poisson 
fluctuations and the contribution from the sampling vari¬ 
ance. The latter is due to the field-to-field variations arising 
from source clustering and is important especially in the case 
of surveys covering small areas. The total fractional variance 
of the differential counts, n, can be written as (|Peebleslll980l ): 


n ~ i n ) \ 2 
(n) / 


1 2 

7—r + Oci 
(n) 
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with 

a 2 cl = -^ J J w(d)<m 1 dQ 2 (2) 

where (n) is the mean count in the flux density bin, 9 is the 
angle between the solid angle elements dQ ,i and dfl 2 , w(9) 
is the angular correlation function and the integrals are over 
the solid angle, f2, covered by the survey. 

The angular correlation function of the fa int sources 
detect able by the SMI camera is not known. iFang et al.l 
(120081) obtained an angular clustering amplitude A ~ 0.001 
at one degree for the galaxies detected at 24 /an by the 
Spitzer Wide-Area Infrared Extragalactic Survey (SWIRE) 
with S' 24 /im > 350 /iJy. The fainter galaxies detected by 
the SMI camera are presumably less luminous on average 
and therefore are unlikely to have a higher clustering am¬ 
plitude. The slope of the 24 /im angular correlation func¬ 
tion is uncertain. Adopting the st andard va lue of —0.8 (i.e. 
w{9 ) ~ O.OOl(0/deg) -0 ’ 8 ) we get rtde Zotti et al. IH): 

a 2 cX = 2.36 x 10 -3 (fi/deg 2 ) -0 ' 4 . (3) 

In the case of a survey over a single FoV = 0.017. This 
is a significant, although minor, contribution to the error 
budget for the proposed deep survey. For the shallow survey, 
the contribution due to the sampling variance is negligible. 

A big plus of the SMI camera is its full coverage of 
the 20-37 /im spectral range with R = 20 re solution. Us- 
ing th e PAH emission template described in iGroves et al] 
J2008I . their eq. (7)) and the normalized parameters of the 
Lorentzian components of the PAH emission band given in 
their Table 2 we have verified that the four PAH bands we 
consider fill the spectral resolution element of the camera 
for the whole redshift range of interest, so that their signal 
is not diluted. In the worst case the fraction of the PAH line 
flux falling within a resolution element varies from ~ 0.8 
(PAH 7.7 and 8.6 /an) to ~ 1 (PAH 6.2 and 11.3 /im). 

Coupling the relationships between the PAH and the 
IR luminosities discussed in Sect.[|]with the redshift depen¬ 
dent IR luminosity functions given by the model we have 
computed the integral counts of galaxies in the four PAH 
lines over the wavelength range covered by the SMI wide 
field camera (see Fig. 0. We find that the proposed shallow 
survey will detect ~ 3.0 x 10 6 galaxies in at least one PAH 
line and ~ 1.9 x 10 6 in at least two lines; for the deep sur¬ 
vey the number of detections are ~ 100 in at least one line 
and ~ 20 in at least two lines. The redshift distributions of 
galaxies detected in 1, 2, 3 and 4 lines are shown in Fig. 0 

Figure [5] shows the minimum SFR (calculated us- 
ing our line/Lm relatio ns and the Air/SFR relation by 
iKennicutt fe Evanj20lJ 'l of the sources detectable (in imag¬ 
ing and in spectroscopy) by the proposed deep survey as a 
function of the redshift. Also shown, for comparison, are 
the SFRs associated to the minimum luminosities repre¬ 
sented in the IR luminosity function s at several redshifts 
determined bv lGruppioni et al.l (hobt ) on the basis of Her- 
schel /PACS and SPIRE surveys. The improvement over 
Herschel is impressive. The deep survey will sample SFRs 
well below those of t he most efficient star formers, estimated 
to be ~ lOOMg/yr dForster Schreiber et al.|[2006l : ICai et al.l 
l2013l l. It will therefore allow a full reconstruction of the dust 
obscured cosmic star formation history up to high redshifts. 
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Figure 5. Comparison between the minimum SFR achieved 
by the proposed deep survey with the SMI wide field camera 
(through spectroscopic detections of PAH lines) and the SFR cor¬ 
responding to the minimum luminosities represented in the IR lu- 
minositv functions determined bvlGruDDioni et al.l J2013I1 on the 
basis of Herschel/PACS and SPIRE surveys, as a function of the 
redshift. 
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H2-5.51/un 

-4.29 

0.39 

PAH6.2/un 2 

-2.20 

0.36 

H26.91/nn 

-3.97 

0.39 

[ArII]6.98/tm 

-3.96 

0.32 

PAH7.7/un 2 

-1.64 

0.36 

PAH8.6/un 2 

-2.16 

0.36 

[ArIII]8.99/i,m 

-4.22 

0.69 

H29.66/mi 2 

-3.96 

0.52 

[SIV] 10.49/un 2 

-3.95 

0.69 

PAH11.3/UTL 1 

-2.29 

0.36 

H 2 12.28/nn 2 

-4.12 

0.54 

HI12.37/xm 

-4.85 

0.30 

[Nell] 12.81/un 1 

-3.11 

0.45 

[ClII]14.38 i um 

-5.44 

0.33 

[Neill] 15.55/un 1 

-3.69 

0.47 

H2l7.03/im 1 

-4.04 

0.46 

[FeII]17.93^m 

-5.18 

0.34 

[SIIIJlSTl/tm 1 

-3.49 

0.48 

[ArIII]21.82//m 

-5.40 

0.70 

[FeIII]22.90/un 

-6.56 

0.33 

[FeII]25.98^im 

-4.29 

0.44 

[SHIJSS^S/im 1 

-3.05 

0.31 

[SiII]34.82 i um 1 

-2.91 

0.28 


1 Taken from iBoimt.n et. aT 1 <2014a|) 
^ Taken from lBonato et al 1 <2014b|) 


Table 1. Mean values of the log of line to IR (8-1000 /^m) con¬ 
tinuum luminosities for star-forming galaxies, (log(L^/Z/iR)), and 
associated dispersions a. For the PAH 1 1.3 and H 2 17.03 fim 
lines, the tabulated values, taken from lBonato et alj d2014alh have 
been computed excluding local ULIRGs, for which the luminosity 
in these lines was found to be uncorrelated with Lir; the mean 
values of log(L^) in these two lines for local ULIRGs, log(L^/Z/Q), 
were found to be 9.01 and 8.07 (with dispersions of 0.28 and 0.34) 
respectively. 
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Figure 6. Comparison between observed (blue squares) and exp ected (red circles) line lu minosities, obta i ned u sing ITERA (see Section[4] 
for details). The [NeVI]7.63 pm data (for AGNs) comes from the lsturm et al] d2Q02f) andlVeilleux et~ al.l d2009l) catalogues. For the other 
3 lines the observed line luminosities (for star-forming galaxies) are from lWillett et al.il d201lh . The green bands show the d=2 <j dispersion 
around the line to continuum luminosity relations derived with ITERA. 
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Figu re 7. Strength o f the 9.7 and 18.0/xm silicate b a nds, as a f unction of the IR luminosity, for star-f orming galaxies. Data from 
IStierwalt et al.1 <|2Q13h . flmanishi et al.1 {2007]), [i manishil {20091 ). and llmanishi. Maiolino. &; Nakagawal (|2010l ). The green bands show the 
2cr intervals around the mean linear relations s\ = — log(Li^) + c (blue lines). 
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4 LINE VERSUS IR LUMINOSITY 

To estimate the counts of galaxy and AGN line detections 
by SMI surveys we coupled the redshift dependent IR (in 
the case of galaxies) or bolometric (in the case of AGNs) lu¬ 
minosity functions of the source populations with relation¬ 
ships between line and IR or bolometric luminosities. We 
have considered the following set of 41 IR lines: 

• 3 coronal region lines: [MgVIII]3.03, [SilX]3.92 and 
[SiVII]6.50pm; 

• 13 AGN fine-structure emission lines: [CaIV]3.21, 
[CaV]4.20, [MgIV]4.49, [ArVI]4.52, [MgV]5.60, [NeVI]7.63, 
[ArV]7.90, [CaV] 11.48, [ArV]13.09, [MgV] 13.50, [NeV] 14.32, 
[NeV]24.31 and [OIV]25.89pm; 

• 19 fine-structure emission lines that can also be pro¬ 
duced in star-formation regions: 

- 10 stellar/HII region lines: [ArII]6.98, [ArIII]8.99, 

[SIV] 10.49, HI12.37, [Nell] 12.81, [GUI] 14.38, [Neill] 15.55, 

[Sill] 18.71, [ArIII]21.82 and [SIII]33.48pm; 


- 4 lines from photodissociation regions: [Fell] 17.93, 
[FeIII]22.90, [FeII]25.98 and [Sill] 34.82 pm; 

- 5 molecular hydrogen lines: H 2 5.51, H 2 6.91, H 2 9.66, 
H 2 12.28 and H 2 17.03 pm; 

• 4 Polycyclic Aromatic Hydrocarbon (PAH) lines at 6.2, 
7.7, 8.6 and 11.3 pm; 

• the 2 emission and absorption silicate bands at 9.7 and 
18.0 pm. 


For the PAH 6.2, PAH 7.7, PAH 8.6, PAH 11.3, H 2 9.66, 
[SIV] 10.49, H 2 12.28, [NeII]12.81, [NeV] 14.32, [Neill] 15.55, 
H 2 17.03, [Sill] 18.71, [NeV]24.31, [OIV]25.89, [SIIIJ33.48 
and [SiI I]34.82 pm lin e s we h ave used the relationships de¬ 
rived bv lBonato et al.l (|2014al lbl) on the basis of observations 
collected from the literature. 

For all the other lines, with either missing or insuffi¬ 
cient data, the line to continuum luminosity relations were 
derived using the IDL Tool for Emission-line Ratio Analysis 
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Figure 8. Strength of t he 9.7 and 18.0 /im sil icate bands as a function of the bolometric luminosity for type 1 (circles) and type 2 
(squares) AGNs from the lGallimore et alj <|2010|) sample. The green or orange bands show the 2cr spreads around the mean strength for 
typel’s (yellow lines) or type2’s (blue lines), respectively. 



Figure 9. Predicted luminosity functions of the AGN lines [NeVI] 7.63, [ArV] 13.09, and [MgV] 13.50 fim at different redshifts at which 
these lines can be detected by the SMI spectrometers. 



Figure 10. Examples of predicted luminosity functions for the total (starburst plus AGN) emission (solid lines) and for the AGN 
component only (dotted lines) of some lines of our sample, at different redshifts at which these lines can be detected by the SMI 
spectrometers. 
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Figure 11. Integral counts in 30 MIR lines of galaxies as a whole (starburst plus AGN components; solid lines) and of AGNs only 
(dotted lines) over the Spec-S and the Spec-L wavelength ranges. The vertical lines correspond to the detection limits (2.0 ■ 10“ 19 W/m 2 
for Spec-S and 8.0 ■ 10 -19 W/m 2 for Spec-L) for 1-h exposure per FoV. 


(ITERA written by Brent Groves. ITERA uses the library 
of published photoionization and shock models for line emis¬ 
sion of astrophysical plasmas produced by the Modelling 
And Prediction in Photolonised Nebulae and Gasdynamical 
Shocks (MAPPINGS III) code. 

Among the options off ered by ITE RA w e have 
chosen, for starbursts, the iDopita et alJ (120061 ) mod¬ 
els and, for AGNs, the dust free isochoric narrow 
line region (NLR) models for typel’s and the dusty 
radiation-pressure dominated NLR models for type2’s 
dGroves, Dopita. fe Sutherland! 120041 '). The chosen models 
are those which provide the best overall fit (minimum 


3 http://home.strw.leidenuniv.nl/ brent/itera.html 


y 2 ) to the observed line ratios of local starbursts in the 
I Bernard-Salas et al. ( 2 0091) cat a logue a nd of AGNs in the 
sample built by Bonatp_eL_aL| ( 2014bT) combining source s 
from the Tommasmet_aL ( 20081 . 20ld ). ISturm et all d2002h . 
and lYeilleux et al.l 1 20091) catalogues. 

As an example we compare, in Fig. [6] the line lu¬ 
minosities as a function of IR luminosities obtained us¬ 
ing ITERA with the o bserved ones (the majority of w hich 
wer e published after the iGroves. Dopita. fe Sutherlandll2004l 
and lDopita et ai1l200d models) for [NeVI]7.63, [ArIII]8.99, 
HI12.37 and [Fell] 17.93 /im. As a further test, we have 
compared the lumi nosities of the f ainter high-ionization 
lines measured by ISpinoglio et al.1 (120051) in the proto¬ 
type Seyfert 2 galaxy NGC 1068 with those obtained 
via ITERA. Adopting for the active nucleus of this ob- 
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Spectral line 

a 

b 

disp (lcr) 

[MgVIII]3.03/nn 

0.73 

-1.65 

0.60 

[CaIV]3.21/im 

0.89 

-2.90 

0.31 

[SiIX]3.92/im 

0.65 

-1.77 

0.88 

[CaV]4.20/rm 

0.90 

-4.00 

0.34 

[MgIV]4.49/un 

0.89 

-3.01 

0.32 

[ArVI]4.52//m 

0.83 

-3.36 

0.34 

H25.51/an 

1.23 

-7.01 

0.40 

[MgV]5.60/nn 

0.91 

-2.94 

0.34 

[SiVII]6.50/im 

0.83 

-3.55 

0.37 

H26.91/-im 

0.80 

-2.40 

0.34 

[ArII]6.98/im 

0.84 

-4.21 

0.64 

[NeVI]7.63/nn 

0.79 

-1.48 

0.42 

[ArV]7.90/an 

0.87 

-3.85 

0.32 

[ArIII]8.99/nn 

0.98 

-4.15 

0.37 

H 2 9.66/un 1 

1.07 

-5.32 

0.34 

[SIVJlOTg/im 1 

0.90 

-2.96 

0.24 

[CaV]11.48/im 

0.90 

-5.12 

0.34 

H 2 12.28/im 1 

0.94 

-3.88 

0.24 

HI12.37/mi 

0.86 

-3.84 

0.34 

[Nell] 12.81/un 1 

0.98 

-4.06 

0.37 

[ArV]13.09/un 

0.87 

-3.85 

0.32 

[Mg V] 13.50/un 

0.91 

-4.01 

0.34 

[NeV]14.32/mi 1 

0.78 

-1.61 

0.39 

[ClII]14.38jam 

0.85 

-4.83 

0.57 

[NelUJlSAS/tm 1 

0.78 

-1.44 

0.31 

11217.03/tm 1 

1.05 

-5.10 

0.42 

[Fell] 17.93/an 

0.84 

-3.80 

0.54 

[3111118.71/im 1 

0.96 

-3.75 

0.31 

[ArIII]21.82/xm 

0.98 

-5.34 

0.36 

[FeIII]22.90/tm 

0.79 

-4.85 

0.60 

[NeV]24.31/mi 1 

0.69 

-0.50 

0.39 

[OIV]25.89/tm 1 

0.70 

-0.04 

0.42 

[FeII]25.98^m 

0.87 

-3.71 

0.55 

[3111133.48/mi 1 

0.62 

0.35 

0.30 

[SiII]34.82p,m 

0.89 

-3.14 

0.52 


1 Taken from |Bonato et al.| d2014b|) 


Table 2. Coefficients of the best-fit linear relations between line 
and AGN bolometric luminosities, log(L^) = a- log(Lboi) + 6, and 
lcr dispersions associated to the relations. 


ject a bolometric luminosity of 3 x 10 11 L@ (iBock et al.l 
I 2000 I) we find, from ITERA, log(L/) in the ranges 
(2cr) [6.13,7.33], [4.81,6.57], [6.88,7.52], [5.83,6.51], [7.16,7.84] 
for [MgVIII]3.03, [SiIX]3.92, [MgIV]4.49, [ArVI]4.52 and 
[MgV]5.60/im, respectively, in reasonably good agreement 
with the corresponding measured values (6.95, 6.60, 6.79, 
7.08 and 7.16, respectively), especially taking into account 
the substantial uncertainty in the estimated bolometric lu¬ 
minosity. There is no indication of systematic over- or under¬ 
estimate of the line luminosities. 

Note that our counts take into account the absorption 
of the fine-structure lines near the strong silicate absorp¬ 
tion features, since such absorption is properly dealt with 
by ITERA. 

The data on starburst galaxies are consistent with a 
direct proportionality between line and IR luminosity. The 
mean line to IR luminosity ratios, (\og(Le/Lni)), and the 
dispersions, a, around them are listed in Tablc[l] In the case 
of AGNs the data are described by linear mean relations 
log(L/) = a ■ log(Lboi) + b. The coefficients of such relations 
and the dispersions around them are listed in Table [2] 


From a practical viewpoint the faintest lines considered 
above will hardly be detected by SPICA instruments. Nev¬ 
ertheless we thought it useful to derive relationships with 
the continuum luminosity for as many IR lines as possible 
to provide a tool to estimate exposure times to detect such 
lines with pointed observations, not necessarily only with 
SPICA instruments. 

For the 9.7 and 18.0 /jm silicate bands we have used 
the observed correlations between the IR luminosity and 
the relative stren gth of the features, defined (see, e.g., 
ISpoon et al.ll2007h as the natural logarithm of the ratio be¬ 
tween the observed flux density at the center of the silicate 
feature, F a b s ,Ay(i+z)i and the local continuum flux density, 

Pcont,A f (1 + z) , 


-Fobs, A f (1+z) 

= ln p -- 

” cont, X f (1+z) 


In 


Fobs, Ay 
Fcont, A f 


(4) 


where A/ is the rest-frame wavelength of the feature (i.e. 9.7 
or 18 pm), while L 0 h S} x f and L con t } x f are the corresponding 
(monochromatic) luminosities at that wavelength. 

To calibrate the relationships between the silicate band 
strength and the IR luminosity for the starb u rst co mpo- 
nent we have used data from Istierwalt et al.l d2013l . only 
for 9.7/nn silicate band), excluding the objects with low 
6.2 fj,m PAH equivalent widths (EQW 6 2/im < 0.27/im) 
whose MIR emission is likely to be substantially con¬ 
taminated b y an AGN, and the starburst dom inated 
gal axies fr om j lmanishi et al.l (i200 7f). Jlmanis hj| (|200flh . and 
llmanishi. Maiolino. fe NakaeawT i 201C)h catalogues (these 
authors actually provide optical depths, tx = —sa). 

As illustrated by Fig. [7] the 9.7 /im strength of starburst 
galaxies shows a clear linear anti-correlation with the log of 
the IR (8-1000/rm) luminosity, consistent with a constant 
c = sa + log Lm with a mean value (c) = 10.92 and dis¬ 
persion of 0.37. The linear correlation coefficient is —0.78, 
corresponding to a correlation significant at the 15.7 a level. 

The very few data on the 18.0 /tm strength do not show 
any signihcant correlation with Lir. However, given the poor 
statistics, the possibility of a correlation cannot be ruled out 
either. If, in analogy to the 9.7 /tm strength, we assume also 
for the 18.0 /tm one a relation of the form s\ = — log Lir + c 
we get (c) = 11.61 with a dispersion of 0.18. If, instead, the 
two quantities are uncorrelated the data give (sa) = —0.57 
with a dispersion of 0.11. In the latter case, the number of 
galaxies detectable in that band by the SPICA SMI spec¬ 
trometers in 1 hour integration per FoV decreases by a factor 
~ 6 . 

For AGNs we have used the iGallimore et al.l (|2f)lf)l i 
sample, neglecting the silicate absorption for typel’s and 
the emission for type2’s. As illustrated by Fig. [8] the sil¬ 
icate strengths of AGNs appear to be uncorrelated with 
the bolometric luminosities. We have therefore adopted 
Gaussian distributions of s\ around mean values, (s\), 
independent of Lboi- We have obtained: (s 9 . 7 ^ m )AGN l = 
0.18, (s9.7^ m )AGN2 = —0.26, (sis.o^m) agn i = 0.08 and 
(si 8 . 0 Mm)AGN 2 = —0.13 with dispersions of 0.14, 0.24, 0.05, 
and 0.12, respectively. 

The line luminosity functions have been computed 
starting from the redshift-dependent IR luminosity func¬ 
tions given by the evolutionary model, including both the 
starburst and the AGN component. To properly take into 
account the dispersion in the relationships between line and 
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Figure 12. Contributions of different AGN populations to the 
SMI Spec-S (solid lines) and Spec-L (dotted lines) integral counts 
in 6 AGN lines ([NeVI]7.63, [ArV]13.09, [MgV] 13.50, [NeV]14.32, 
[NeV]24.31 and [OTVI25.S9//rri). A: AGNs associated to late- 
type galaxies plus optically selected AGNs; B: AGNs associ¬ 
ated to unlensed proto-spheroids; C: AGNs associated to strongly 
gravitationally lensed proto-spheroids. The vertical lines corre¬ 
spond to the detection limits (2.0 • 10“ 19 W/m 2 for Spec-S and 
8.0 ■ 10 -19 W/m 2 for Spec-L) for 1-h exposure per FoV. 
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Figure 13. Predicted redshift distributions of AGNs detected in 
6 lines ([NeVI]7.63, [ArV]13.09, [MgV]13.50, [NeV] 14.32, [NeV]- 
24.31 and [OIV]25.89 pm) detected by SMI Spec-S (solid red lines) 
and Spec-L (dashed blue lines), compared, for three of them 
([ NeV]14.32, [NeVl24.3 1 and [OIVJ25.89 pm), with those obtained 
in lBonato et al.l (|2(I 1 111 ) for SPICA/SAFARI (dotted black lines), 
for an 1 h exposure per FoV. 


continuum luminosi ties we have used th e Monte Carlo ap¬ 
proach described in lBonato et al.l d2014bl ). Examples of line 
luminosity functions at various redshifts are shown in Figs. [9] 
and EH] 


5 LINE LUMINOSITY FUNCTIONS AND 
NUMBER COUNTS 

Our predictions for the integral counts in both the Spec-S 
and the Spec-L channels, for 30 lines are shown in Fig. 1111 In 
1 li integration per FoV, the Spec-S will detect ~ 70 AGN 
lines per square degree, primarily the [NeV]24.31 pm line 
(~ 40% of detections); the Spec-L will detect mostly the 
[OIV]25.89pm line (over 90% of the ~ 70 AGN line detec¬ 
tions). The ~ 140 lines are produced by ~ 110 individual 
AGNs. Only a handful of them are optically selected. Thus 
this survey will be very efficient at selecting obscured AGNs, 
complementing optical surveys. 


Figure fl2l shows the contributions of AGNs associated 
with different galaxy populations to the SMI Spec-S and 
Spec-L counts in six AGN lines. The redshift distributions of 
AGNs detected in each of these lines in 1 li integration/FoV 
are displayed in Fig. [TS] For three of the lines ([NeV] 14.32, 
[NeV]24.31 and [OIV]25.89 pm) we also show, f or compar- 
ison, the redshift distributions obtained by iBonato et al.l 
(l2014bf) for a SPICA/SAFARI survey (again for a lh ex¬ 
posure per FoV). The two instruments cover nicely comple¬ 
mentary redshift intervals. 

Figure[TT]illustrates the redshift distributions of sources 
for which the SMI spectrometers will detect at least 1, 2, 3 or 
4 lines in 1 h integration per FoV over an area of 1 deg 2 . The 
corresponding numbers of detected lines are about 53,000, 
16,000, 3,100, 360, 55 and 8, respectively. Sources detected 
in at least one line include ~ 200 strongly lensed galaxies at 
z > 1. 

Note that the 5 a detection of only one line does not 
necessarily imply that the redshift determination is prob¬ 
lematic. Many of the objects detected in only one line can 
show a suite of weak but measurable features (other lines, 
absorptions, PAH bumps). The global pattern can then allow 
the determination of reliable redshifts even when individual 
features are only significant at the 2-3 a level. 

Some multiple line detections can be due to different 
galaxies seen by chance within the same resolution element. 
We have estimated this confusion effect adopting a typical 
FWHM of 2" and assuming a random galaxy distribution 
(negligible clustering effects). The fractions of two line de¬ 
tections due to confusion by the SMI camera and by the 
SMI spectrometers are shown as a function of the integra¬ 
tion time per FoV in Fig. 1151 The confused fraction is always 
small (in particular in the case of SMI spectrometers). For 
example, with an integration time of 1 hr two line detec¬ 
tions due to confusion are ~ 1.5% (for the camera survey) 
and ~ 0.2% (for the spectrometers). The number of confu¬ 
sion cases grows almost linearly with the integration time. 
Therefore the confused fraction grows as we go to fainter 
fluxes, where the number counts are flatter, but it is still 
only ~ 3.6% (camera) and ~ 0.7% (spectrometers) for an 
integration time of 10 h per FoV. 

The difference between the AGN and the galaxy SEDs 
in the SMI range implies that the equivalent widths (EQWs) 
of the brightest spectral lines excited by star formation are 
useful indicators of the AGN contribution. This is illustrated 
by Fig. [16] which shows the variation of the EQWs of the 
most prominent star formation lines (log(EQW) > —1.5) 
with the fractional AGN contribution to the total (starburst 
plus AGN) Lir( 8-1000/im). The PAH lines are particularly 
effective for this purpose. 


6 OBSERVING STRATEGY 

As illustrated by Fig. 1111 the integral counts for both SMI 
spectrometers have a slope flatter than 2 at and below the 
detection limit for 1 h integration/FoV for the majority of 
the lines. Counts of PAH lines with the SMI camera show 
a similar behaviour (Fig. [3j. This means that the number 
of detections for a fixed observing time generally increases 
more by extending the survey area than by going deeper, 
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Figure 14. Predicted redshift distributions of galaxies (starburst 
plus AGN components) detectable in one (black histogram), two 
(cyan), three (red) and four (blue) spectral lines, by a SPICA 
SMI spectrometer (Spec-S plus Spec-L) survey covering 1 deg 2 in 
1 h integration/FoV. 


Figure 15. Detection of multiple spectral lines produced by 
galaxies at different redshifts falling within the same resolution 
element: the fraction of cases of confused lines is shown as a func¬ 
tion of the integration time per FoV, for a survey with the SMI 
camera (solid blue line) and with the SMI spectrometers (dotted 
red line). 


similarly to what found bv lBonato et alJ (l2014al lbh for blind 
spectroscopic surveys with SPICA/SAFARI. 

As mentioned in Sect. [5] we expect that a survey of 
1 deg 2 with 111 integration/FoV will detect ~ 110 AGNs. 
Therefore, to investigate the AGN evolution with suffi¬ 
cient statistics we need a much wider-area. Also, the blind 
SMI spectroscopic survey may be usefully complemented by 
follow-up observations of bright high -2 galaxies already dis¬ 
covered at (sub-)mm wavelengths over much larger areas. 
Figure [T7] shows the SMI spectrometer exposure time per 
FoV needed to achieve a 5cr detection of typical AGN lines 
at 2 = 1 and 2 = 1.5 as a function of the bolometric lumi¬ 
nosity. We can see, for example, that the [NeV] 14.32 pm line 
can be detected in 10 h from an AGN with (real or apparent, 
i.e. boosted by strong gravitational lensing) bolometric lu¬ 
minosity of ~ 10 13 Lq at 2 = 1 and ~ 4 x 10 13 Lq at 2 = 1.5. 
When the South Pole Telescope (SPT) and the Herschel sur¬ 
vey data will be fully available we will have samples of many 
hundreds of galaxies with either intrinsic or apparent IR lu- 
minosities l arger than 10 13 Lq llNegrello et al]l 201i l 120141: 

I Vieira et al.ll2013l ). As explained in lBonato et alJ ( 2014bf) . 
pointed observations of those sources can allow us to inves¬ 
tigate early phases of the galaxy/AGN co-evolution. 


7 CONCLUSIONS 

We have worked out predictions for surveys with the SMI 
wide field camera and with the spectrometers. 

The combination of a shallow and of a deep survey with 
the camera, requiring a total observing time of ~ 1046 h, 
will allow an accurate definition of MIR source counts of 
both galaxies and AGNs over about five decades in flux den¬ 
sity, down to ~ 3.0 /rJy, i.e. more than one order of magni¬ 
tude fainter than the deepest Spitzer surveys at 24 pm. This 
amounts to resolving almost entirely the MIR extragalactic 
background. The spectral resolution of the camera is opti¬ 
mally suited to detect PAH lines, yielding redshift measure¬ 
ments. The redshift information will allow us to derive the 
SFR function down to SFRs hundreds of times lower than 



Figure 16. Equivalent widths (EQWs) of the brightest spectral 
lines in our sample for a galaxy with a fixed total (star forming 
plus AGN) IR luminosity of 10 13 Lq and varying fractions of AGN 
IR luminosity. 


was possible using Herschel surveys and well below the SFRs 
of typical star-forming galaxies. The cosmic dust obscured 
star formation history will then be accurately determined at 
least up to 2 ~ 4. 

On the spectroscopic side we have considered 41 MIR 
lines, 6 of which are predominantly excited by AGN activity 
while the others are primarily associated with star forma¬ 
tion. Relationships between the line luminosity and the IR 
(for the starburst component) and/or the bolometric lumi¬ 
nosity (for the AGN component) are presented. Several of 
them were derived in previous papers, but many are new. 

Using these relationships we computed the expected 
number counts for the 30 brightest lines. We found that 
the SMI spectrometers will detect, with an integration time 
of lh/FoV, about 52,000 galaxies per square degree in at 
least one line and about 16,000 in at least two lines. About 
200 of galaxies detected in at least one line will be strongly 
lensed. 

The number of expected AGN detections is far lower. 
For the same integration time (lh/FoV) we expect to de¬ 
tect, per square degree, ~ 140 AGN lines, from ~ 110 AGNs. 
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Figure 17. SPICA SMI Spec-S (solid lines) and Spec-L (dotted 
lines) exposure time per FoV required for a on line detection of 
3 typical AGN lines ([ArV]13.09, [MgV]13.50 and [NeV] 14.32/un) 
as a function of the AGN bolometric luminosity for z = 1.0 (left) 
and z = 1.5 (right). 
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Thus a larger area is necessary to investigate the AGN evo¬ 
lution with good statistics. 

Given the low surface density of AGNs detectable by 
the SMI spectrometers, an efficient way to investigate early 
phases of the galaxy/AGN co-evolution are pointed observa¬ 
tions of the brightest galaxies detected by large area surveys 
such as those by Herschel and by the SPT. 
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